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ABSTRACT 

We use the IRAM HERACLES survey to study CO emission from 33 nearby spiral galaxies down to 
very low intensities. Using 21-cm line atomic hydrogen (Hi) data, mostly from THINGS, we predict 
the local mean CO velocity based on the mean H I velocity. By re-normalizing the CO velocity axis 
so that zero corresponds to the local mean H I velocity we are able to stack spectra coherently over 
large regions. This enables us to measure CO intensities with high significance as low as Iqo ~ 0.3 
K km s -1 (Sh2 ~ 1 Mq pc~ 2 ), an improvement of about one order of magnitude over previous 
studies. We detect CO out to galactocentric radii r ga i ~ r25 and find the CO radial profile to follow 
a remarkably uniform exponential decline with scale length of ^0.2 T25. Here we focus on stacking 
as a function of radius, comparing our sensitive CO profiles to matched profiles of Hi, Ha, FUV, 
and IR emission at 24 /im and 70 /xm. We observe a tight, roughly linear relationship between 
CO and IR intensity that does not show any notable break between regions that are dominated by 
molecular gas (£h2 > Shi) and those dominated by atomic gas (Sh2 < Shi). We use combinations of 
FUV+24/xm and Ha+24/xm to estimate the recent star formation rate (SFR) surface density, Ssfr, 
and find approximately linear relations between Esfr and Eh2 • We interpret this as evidence for stars 
forming in molecular gas with little dependence on the local total gas surface density. While galaxies 
display small internal variations in the SFR-to-Fb ratio, we do observe systematic galaxy-to-galaxy 
variations. These galaxy-to-galaxy variations dominate the scatter in relationships between CO and 
SFR tracers measured at large scales. The variations have the sense that less massive galaxies exhibit 
larger ratios of SFR-to-CO than massive galaxies. Unlike the SFR-to-CO ratio, the balance between 
atomic and molecular gas depends strongly on the total gas surface density and galactocentric radius. 
It must also depend on additional parameters. Our results reinforce and extend to lower surface 
densities a picture in which star formation in galaxies is separable into two processes: the assembly 
of star-forming molecular clouds and the formation of stars from H2. The interplay between these 
processes yields a total gas-SFR relation with a changing slope, which has previously been observed 
and identified as a star formation threshold. 

Subject headings: galaxies: evolution — galaxies: ISM — radio lines: galaxies — stars: formation 



1. INTRODUCTION 

Stars form out of molecular (H2) gas and many re- 
cent observations of nearby galaxies have revealed a 
strong correlation between the surface density of molec- 
ular gas, Eh2, and the star formation rate (SFR) sur- 
face density, E SF r fWong fc BlitzH2002t iKennicutt et al.l 
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2008; 



Bigiel et aJJlpOOl iLerovetaLl [20081: IWilson et all 
Blanc et al.l 2009; see also the recent review by 



Bigi el et alJl20TTF . These stud es show a correlation over 



several orders of magnitude, but mostly for regions where 
H2 makes up the majority of the neutral gas, Sh2 ^ Em- 
The lack of a clear correlation between atomic gas (H 1) 
surface density, S hi, and Ssfr inside galaxy disks (e.g., 
iBigiel et al.ll2008i ) offers circumstantial evidence that star 
formation remains coupled to the molecular, rather than 
total (H1+H2), gas surface density even where Hi makes 
up most of the interstellar medium (ISM). However, the 
exact relationship between Ssfr and Sh2 in the Hi- 
dominated parts of galaxies (Sh2 ^ Shi) remains largely 
unexplored. 

In this paper we use new, sensitive, wi de-field CO 
maps from the IRAMlEI HERACLES survey (|Lerov et al.l 
120091) to measure correlations between molecular gas and 
SFR tracers over a large dynamic range. By employing 
stacking techniques based on Hi priors we extend our 
observations from H2 -dominated galaxy centers to the 
outer parts of galaxies where the H2 surface density is 

12 IRAM is supported by CNRS/INSU (France), the MPG (Ger- 
many) and the IGN (Spain). 
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much lower than the Hi surface density, Ejj2 Shi- 

Deep CO measurements allow us to test if a single 
"star formation law" applies in both t he H2- and Hi - 
dominated parts of galaxies. Following iSchmidtl (|1959f ) . 
astronomers have investigated scaling relations linking 
gas and star formation for decades. Such relations only 
approximate the complex physical processes involved in 
star formation but provide useful constraints on theoret- 
ical model s and impor tant input to simulations. After 
IKennicuttl ffl989l . fl998h . power laws linking surface den- 
sities of gas and SFR are the most common formulation. 
However, the choice of which gas surface density to use — 
total or molecular gas — remains controversial, as does 
the extension of any measured molecular relation to low 
surface densities. The underlying question is what limits 
star formation in low column density regions, the for- 
mation of molecular gas or the efficiency at which the 
available molecular gas is converted into stars? Sensitive 
observations of molecular gas down to low surface densi- 
ties in a large sample of galaxies are needed to address 
these questions. 

Our CO measurements also allow us to investigate the 
distribution of molecular gas out to large radii. A char- 
acteris tic exponential decline has been observed several 
times (|Young et al.i ri995; Re gan et al.ll2001t iLerov et al.l 
2009) , but it is not known if this decline becomes sharper 
at one point, for example c orresponding to claimed sta r 
formation thresholds (e.g., iMartin fc Kennicuttl [2001). 
We also test how variations in the H2-to-Hi ratio 
extend to low surface densities. This quantity is a 
strong a nd systematic function of environment in nearby 
galaxies flWong fc Blitz! [2001 iBlitz fc Rosolowsk^ [2006l: 
ILerov et al.ll200& Iffitschfeld et al.ll2009D but it has been 
difficult to extend the observed correlations to low sur- 
face densities. 

In Section [2] we describe our sample and data. In Sec- 
tion [3] we present the method that we use to extract sen- 
sitive CO measurements. In Section 0] we present radial 
profiles of Hi, CO, IR, FUV, and Ha and use these to 
relate CO, Hi, and tracers of recent star formation. In 
Section [5] we summarize our results. 

2. SAMPLE & DATA 

We study 33 nearby, star-forming disk galaxies, the 
set of HERACLES targets for which we could collect the 
necessary Hi, IR, FUV, and Ha data. This is mainly 
the overlap of several surveys: HERACLES (IRAM 30m 
CO . ILerov et al.l[2009h . THINGS (VLA Hi.lWalter et al l 
20081). SINGSor LVL CffpifeerlR. lKennicutt et aTT2 003a: 
Dale et al.ll2009h and 1 the GALEX NGS (GALEX FUV, 
Gil de Paz et al.l 12007ft . We supplement these with a 
combination of archival and new Hi data and archival 
GALEX data. We exclude low mass, low metallicity 
galaxies with only upper limits on CO emission and 
nearly edge-on galaxies. 

Table Q] lists our sample along with adopted mor- 
phology, distance, inclination, position angle, optical 
radius, and metalli city. These values are taken from 
iWalter et all (|2008f ) if possible and from LEDA and 
NED in other case s. We quote oxygen abundances 
(metallicities) from iMoustakas et al.l (|2010t Table 9), 
averaging the metallicities derived from a theoretical 
calibration (their KK04 values) and an empirical cal- 
ibration (their PT05 values). For galaxies without a 
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NGC337 


SBd 


24.7 


51 


90 


1.48 


8.51 


NGC628 a ' b 


Sc 


7.3 
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20 


4.92 


8.68 


NGC925 a > b 


SBcd 


9.2 


66 


287 


5.32 


8.52 


NGC2146 


SBab 


12.8 


54 


123 


2.69 


8.68 d 


NGC2403 a 


SBc 


3.2 


63 


121 


7.87 


8.57 


NGC2798 a 


SBa 


24.7 


85 


152 


1.20 


8.69 


NGC2841 a - b 


Sb 


14.1 


71 


153 


3.45 


8.87 


NGC2903 a - b 


SBd 


8.9 


65 


204 


5.92 


8.90 


NGC2976 a ' b 


Sc 


3.(5 


65 


335 


3.60 


8.67 


NGC3049 


SBab 


8.9 


58 


28 


1.04 


8.82 


NGC3077 a 


Sd 


3.8 
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15 


2.70 


8.64 
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SBc 


13.8 


72 
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NGC3351 a,b 


SBb 


10.1 


11 


192 


3.60 


8.90 


NGC3521 a - b 


SBbc 


10.7 


73 


340 


4.16 


8.70 


NGC3627 


SBb 


9.3 


62 


173 


5.14 


8.67 


"\T /~i t~1 OA'IO 

1NCj03938 


Sc 


12.2 


14 


15 


1.77 


8.74 


NGC4214 a > b 


In 


2.9 


11 


65 


3.40 


8.25 


NGC4254 


Sc 


20.0 


32 


55 


2.51 


8.79 


NGC4321 


SBbc 


14.3 


30 


153 


3.01 


8.83 


IN Vj040v50 
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SBcd 


11.6 


65 


328 


5.24 


8.55 


NGC4569 


SBab 


20.0 


66 


23 


4.56 


8.92 


NGC4579 


SBb 


20.6 


39 


100 


2.51 


8.88 


NGC4625 


SBmp 


9.5 


47 


330 


0.69 


8.70 


NGC4725 


SBab 


9.3 


54 


36 


4.89 


8.73 


NGC4736 a - b 


Sab 


4.7 


41 


296 


3.87 


8.66 


NGC5055 a ' b 


Sbc 


10.1 


59 


102 


5.93 


8.77 


NGC5194 a 


SBc 


8.0 


20 


172 


3.85 


8.86 


NGC5457 a 


SBcd 


7.4 


18 


39 


11.99 


8.46° 


NGC5713 


Scd 


26.5 


48 


11 


1.23 


8.64 


NGC6946 a ' b 


SBc 


5.9 


33 


243 


5.70 


8.72 


NGC7331 a - b 


Scd 


14.7 


76 


168 


4.59 


8.68 



a Targets of THINGS survey l|Walter et al.ll2008h 

Targets in first HERAC LES survey paper IILerov et al.ll2009T) 
c Oxygen abundance from Moustakas et al. (2010) 
d Oxygen abundence from Engclbracht et al. (2008) 
e Oxygen abundance from lKennicutt et a l. (2003b) 



IMoustakas et al.l (|2010l ) metallicity, we adopt a metal- 
licity equal to the average of their B-band luminosity- 
metallicity rel ations. For NGC 2146 we quote the metal- 
licity given bv lEngelbracht et al.l (pOOl ). For NGC 5457 
(M 101) we take a constant metallicit y defined by the 
value a t 0.4 r 2 5 from the gradient fit bv lKennicutt et all 
([200311 . 

We trace molecular hydrogen (H 2 ) using CO (2 — > 1) 
line emission observed w ith the IRAM 30m as part of 
the HERACLES survey (|Lerov et al.l l200^h . They de- 
scribe in detail the observations and reduction for the 
subset of galaxies observed until Summer 2008. The re- 
maining targets were observed and reduced in the same 
way. The final data cubes have an angular resolution 
(FWHM) of 13" and a spectral resolution (channel sep- 
aration) of 2.6 km s — 

Our measurements of atomic h ydrogen (Hi) com e 
mostly from the THINGS survey ((Walter et aLll2008D . 
which used the Very Large ArrayEEl (VLA) to observe the 
21-cm hydrogen line in 34 nearby galaxies. The observ- 
ing and reduction strategies are described therein. The 

13 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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final data cubes have an angular resolution of —11" (us- 
ing natural weighting) and a spectral resolution of 2.6 or 
5.2 km s" 1 . THINGS is sensitive to S H i ~ 0.5 M© pc~ 2 
on scales of 30". Using azimuthal averaging, we reach 
even better sensitivities at large radii. As a result, the 
H I sensitivity never limits our analysis. 

The Hi data for NGC 337, 2146, 2798, 3049, 3938, 
4254, 4321, 4536, 4579, 4625, 4725, and 5713 are a com- 
bination of new and archival VLA data (the new data 
are from VLA programs AL731 and AL735). These have 
been reduced and imaged using the Common Astronomy 
Software Applications (CASA) following a similar pro- 
tocol than the THINGS reduction. These supplemental 
Hi cubes include only data from the VLA's C and D 
configurations; THINGS also includes B configuration 
data. For NGC 4559 we take Hi data observed with 
the Westerbork Syn thesis Radio Telescope (WSRT) by 
Ivan der Huistl (|2002f ). The beam sizes (FWHM) of the 
supplemental H I are 15" — 25" and the velocity resolution 
is 2.5 — 20 km s _1 , usually 10 km s _1 . 

We derive Hi surface densities from 21-cm line inten- 
sities and H2 surface densities from CO(2 — > 1) line in- 
tensities following 

S H i = 0.02 Im x cos i (1) 

Eh2=6.25 Ico x cos i (2) 

where Shi and Em have units of M pc -2 and /hi 
and Ico ar e measured in K km s . The mass sur- 
face densities are projected to face-on values and in- 
clude a factor of 1.36 to account for heavy elements. 
For Equation ([2]) we have assumed a CO line ratio of 
J C o(2 -> l)/Jco(l -> 0) = 0.7 and a CO(l -> 0)-to-H 2 
con version factor Xga = 2.0 x 10 20 cm~ 2 (K km s -1 ) -1 
(see lLerov et al.ll2009L and references therein). 

We use broadband infrared (IR) photometry at 24 /jm 
and 70 /im obtained by t he Spitzer leg acy surveys SING S 
(jKennicutt et all l2003bD and LVL pale et all [2009h . 
Spitzer has angular resolution —6" at 24 /im and 18" 
at 70 /im. The sensitivity of these data is sufficient to 
obtain high signal-to-noise measurements when averag- 
ing in radial rings (see below) throughout the area that 
we study (r gal < 1.2 r 2 s)- 

The GALEX NGS (|Gil de Paz et al.|[200l imaged far- 
and near-ultraviolet (FUV and NUV) emission for most 
of our targets. The FUV band covers 1350 - 1750 A 
with an angular resolution —4.5". We use these images to 
trace unobscured emission from young stars. For galaxies 
not covered by the NGS, we searched the NASA Multi- 
mission Archive at STScI and used the FUV image with 
the longest exposure time. These data also have suffi- 
cient sensitivity to determine FUV intensities with high 
signal-to-noise throughout the star-forming disk. 

We draw Ha data from the SINGS and LVL surveys, 
complemented by literature data for NGC 2903, 4214, 
4569, 4736, and 5457. For the literature and several 
problematic SINGS targets, we pin the total Ha+fNill 
flux to published values, usually those of lKennicutt et all 
(200|). Leroy et al. (2011, in prep.) describe the 
processing of the maps, which involves subtracting a 
smooth background, maski ng foreground stars following 
iMuhoz -Mateo s et al.l (|2009f) , correcting for [N 11] contam- 



ination following iKennicutt et al.l ([2008), and correct- 
ing for Galactic extinction. The Ha maps become un- 
certain, and likely biased low due to background sub- 
traction, below intensities equivalent to Esfr of a few 
times 10 -4 M© yr _1 kpc -2 , values typically crossed in- 
side the radial range studied here. This, and the de- 
clining Ha-to-FUV flux ratios which are observed as 
Galex UV disks extending far beyond t he Ha emission 
(jThilker et alj|2007t iMeurer et atll2009H . limit the util- 
ity of the Ha maps to trace star formation in the low 
brightness regions of outer galaxy disks. 

After we examine correlations among observables, 
we will estimate the star formation rate surface 
density, EsfRi from comb inations of Ha with 
24 um (IKennicutt et all 120071) and FUV with 24 /im 
dLerov et al.ll2008r ) . We adopt the Ha+24/im calibration 
bv ICalzetti et al.1 (120071) and the FUV+24/im combina- 
tion from lLerov et al.l (|2008f ). 

EsFR(Ha+24) = 

(2.9 x 1CT 2 I Ha + 2.5 x 1CT 3 h^m) x cos i (3) 

EsFR(FUV+24) = 

(8.1 x 1(T 2 /fuv + 3.2 x 10~ 3 I 2 ^ m ) x cos i (4) 

where Ssfr has units of M© yr -1 kpc~ 2 and Ha, FUV, 
and 24 /im intensities are all in MJy sr _1 . Both Esfr 
calibrations combine a tracer of the unobscured star for- 
mation with infrared (24 /im) emission, which is intended 
to trace young starlight rep rocessed by dust. Ha traces 
O stars with ages <5 Myrs (Ke nnicutt et al.ll2009f) with 
sensitivity out to -10 Mvr (|Vacca et al.1 11996T K FUV 
traces O and B stars of typi cal age 20 — 30 M yrs with 
sensitivity out to ~ 100 Myr (jSalim et al.ll2007t ). 

Our use of both Ha+24/im and FUV+24/im emission 
gives some test of sensitivity to our choice of star forma- 
tion rate tracer. Several other concerns are worth men- 
tioning. Because dust properties and the stellar popula- 
tions heating the dust somewhat differ between Hll re- 
gions and large (kpc) regions in galaxies the appropriate 
weighting of the 24 /im emission to correct for extinc- 
tion may be a function of scale and environment. The 
dust-to-gas ratio, dust size-distribution, ISM geometry, 
and recent sta r formation his t ory m ay also play impor- 
tant roles. The ICalzetti et al~l (j2007l ) calculation remains 
state of the art, but there is no definitive consensus about 
the correct calibration to use outside bright regions that 
they study. The reliability of the Ha imaging at low 
surface brightness also represents a concern. Ground 
based, narrowband Ha imaging is challenging and the 
Ha+24/tm tracer must be considered unreliable where 
E SF r(hq) ~ 5 x 10 -4 M© yr~ x kpc~ 2 . The Galex and 
Spitzer maps are better behaved at low surface bright- 
ness. 

Because of these data quality considerations and our 
focus on regions with low surface brightness, we em- 
phasize comparisons to FUV+24/im. SFR(FUV+24) 
and SFR(Ha+24) give comparable results with scatter 
of only 0.1 dex (-25%) down to £ S fr - 5 x 1CT 4 
Mq yr -1 kpc~ 2 in an azimuthally averaged ring. Our 
main method to address these other systematic concerns 
is to emphasize the observed scaling relations in the first 
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Before V. Field Alignment After V. Field Alignment 

0.02 I ''I'M 0.021 




-50 50 100 

<v HI > - <v co > [km s _1 ] 

Fig. 1. — The difference between local mean velocities of Hi and 
CO for lines of sight with galactocentric radius smaller than 0.5 rih- 
The close correspondence motivates our use of the H I mean velocity 
to predict the CO velocities in low surface brightness regions. 

part of the paper. Fo r a more thorough dis cussion of 
hybrid SFR tracers see iKennicutt et all (|2009() and for a 
discussion of their application to gas-SFR comparisons, 
we refer the reader to Leroy et al. (2011, in prep.). 

3. METHODOLOGY 

Our goal is to recover low brightness CO emission from 
the outer parts of galaxies. CO is very faint in these 
regions and individual spectra have low signal-to-noise 
ratios (SNR), requiring us to average many spectra to 
achieve a detection. Because the velocity of CO emission 
varies with position, simply averaging spectra spreads 
the emission across many velocity channels with low SNR 
in each channel. In principle, this could still yield a high 
SNR measurement. In practice we wish to maximize 
SNR by considering only the part of each spectrum likely 
to contain emission. We must also contend with sys- 
tematic effects of weather, receiver instabilities, and dish 
imperfections. These all induce frequency-dependent be- 
havior ( "baselines problems" ) that make a clear detection 
of an emission line an important step in a robust analysis. 

With these issues in mind, we use the following tech- 
nique to average CO spectra across large parts of a 
galaxy. First, we estimate the local mean velocity of 
CO emission from the Hi data. Using this mean veloc- 
ity, we redefine the velocity axis of each CO spectrum 
so that the local mean velocity is now zero. We average 
these shifted CO spectra from across our target region. 
In the averaged spectra we expect the CO line to emerge 
at zero velocity with good SNR. The baseline problems 
described above will not average coherently, allowing a 
straightforward identification of the line. 

This approach hinges on the assumption that the Hi 
mean velocity is a good proxy for the mean velocity of the 
molecular gas. Figure Q] shows that this assumption holds 
where we detect CO over individual lines of sight, with 
median vqo — vm of —0.22 km s" 1 and a la dispersion of 
7.0 km s _1 . We expect a similar correspondence in the 
outer, CO-faint parts of galaxies. 

Leveraging Hi to detect CO at large radii works be- 




1 ti , 4 

*v* 4 * * *%*** 



200 400 600 800 
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-400 -200 200 
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Fig. 2. — Average CO spectrum across tilted ring spanning from 
0.7 - 0.8 r 2 5 in NGC 5055. The left panel shows the result of a 
simple average of all spectra. The right panel shows the average 
after each spectrum is shifted so that v = corresponds to the 
local mean H I velocity. 

cause the H I surface densities are essentially constant out 
to large radii, making Hi easily detected across galaxy 
disks. CO emission, on the other hand, tends to be bright 
in galaxy centers but declines rapidly with increasing 
galactocentric radius. 

3.1. Stacking of CO Spectra 

Predicting the velocity of the CO line from the Hi 
data allows us to increase the SNR when measuring the 
integrated CO intensity. The HERACLES bandpass is 
~1000 km s _1 and a typical CO line width at large galac- 
tocentric radii is ^25 km s^ 1 . The Hi allows us to re- 
strict our integration to just the relevant part of the spec- 
trum which represents a substantial gain in sensitivity. 

Just as important as the increase in SNR, the shifted 
and stacked spectra allow us to verify that faint emis- 
sion is actually an astronomical signal. Even for faint 
CO emission the stacking technique has the potential 
to reveal a spectral line. Low-level variations due to 
weather, receiver instabilities, and other systematic ef- 
fects in the telescope will not create such an effect. Even 
stray pickup of astronomical emission due to surface im- 
perfections (i.e., error beam effects) will emerge at a low 
level offset from zero velocity due to galaxy rotation. 

Figure [5] demonstrates this approach. We plot the av- 
eraged CO spectrum of NGC 5055 inside a tilted ring 
spanning from 0.7 — 0.8 r 2 5- In the left panel the spec- 
tra were averaged as they were observed, whereas in the 
right panel we first shifted by the local H I mean velocity 
and then averaged. Both spectra contain the same inte- 
grated intensity, however, only the appearance of a clear 
line feature in the right spectrum at the expected velocity 
strongly indicates that the signal is not due to baseline 
features but cannot be anything but CO emission. 

3.2. Fitting the CO Line 

To extract CO line emission, we perform an automated 
line fit to each stacked spectrum. This approach picks 
out spectral line emission rather than baseline structure 
and does not require us to define an integration window 
beforehand. 

In most regions, the line can be well-approximated by 
a Gaussian profile with FWHM of ^15-40 km s _1 (Fig- 
ure [3] upper right panel). However, in the central regions 
of some galaxies the line can be very broad with a flat- 
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Fig. 3. — Examples of stacked CO spectra with different line 
shapes. The upper row shows spectra with high quality fits and 
the lower row shows spectra with low quality fits. The horizontal 
dotted lines mark the lcr rms noise of the stacked spectrum. The 
upper left panel shows a broad line approximated by a double- 
horn profile. The upper right panel shows a narrow line fitted by 
a Gaussian. The lower left panel shows a marginal quality fit. 
The lower right panel shows a low quality fit together with the 
associated upper limit (dashed line). 

tened or double-horned peak (Figure [3] upper left panel). 
These profiles often coincide with central enhancements 
like bars or molecular rings and are poorly parametrized 
by a single Gaussian. Instead we fit a double-horn pro- 
file, a Gaussian scaled by a sy mmetric second- order poly- 
nomial (functional form from lSainton gc 2007). Based on 
"by eye" inspection, the asymmetry in these profiles is 
small enough that a symmetric function is sufficient. 

We derive the best fit profile via a non-linear least 
squares fi10. We constrain the fit parameters so that the 
center of the profile lies within ±50 km s _1 of zero ve- 
locity after shifting, the FWHM is larger than 15 km s _1 
(to avoid the fit latching onto individual channels), and 
the amplitude is positive. We always carry out a Gaus- 
sian fit first and in those cases where the FWHM exceeds 
60 km s _1 we switch to a double-horn profile. We verify 
by eye that this yields sensible results. 

The integral of the fitted profile gives us the integrated 
CO line intensity. We derive the uncertainty in this quan- 
tity from the noise, estimated from the signal-free part 
of the spectrum, and the width of the profile. It proved 
useful to define a quality scale for the fit. The quality 
is "high" where the peak intensity is larger than 5cr and 
its integrated intensity is larger than 10 times its uncer- 
tainty. In cases where the peak intensity is less than 3cr 
or the integrated intensity is less than 5 times its un- 
certainty, we do not trust the fit and instead determine 
an upper limit. We label cases that fall between these 
regimes as "marginal" . Figure [3] shows examples of our 

14 We use the 1DL procedure MPFIT.PR0 from Craig Markwardt 
which performs a Levenberg-Markwardt non-linear least squares 
minimization and is based on the MINPACK-1 LMDIF.F algorithm 
from More et al. (1978) 



quality measures and line profile fits. 

We derive our upper limits integrating over a Gaussian 
line profile with FWHM set to 18 km s _1 , the typical 
FWHM found for high SNR spectra at r > 0.5 r 25 , and 
fixed amplitude of 3cr. 

3.3. Stacking as a Function of Radius 

We present our stacking technique applied to radial 
bins. In principle this method allows us to stack spec- 
tra across any region. For example, we could define re- 
gions by total gas column, infrared intensity, or features 
such as spiral arms and bars. In practice, radius makes 
an excellent ordinate. We wish to study the underlying 
relationship between CO, Hi, IR, FUV, and Ha inten- 
sity. Stacking with one of these quantities as the ordinate 
would require carefully modeling the biases involved to 
measure the underlying relationships. Galactocentric ra- 
dius is a well-determined, independent quantity that is 
also highly covariant with these other intensities. This 
yields a dataset with large dynamic range that is easy to 
interpret. 

Therefore we focus our analysis on data stacked in bins 
of galactocentric radius. We average over tilted rings 15" 
wide, comparable to the angular resolution of our data. 
This width corresponds to ^220 pc for our nearest tar- 
gets (3 Mpc) and ^1800 pc for our most distant targets 
(25 Mpc). We construct the rings assuming that each 
galaxy is a thin disk with the inclination and position 
angle given in Table [TJ To measure CO with highest 
sensitivity we construct stacked CO spectra using the 
procedure described above and fit those to determine 
the integrated CO line intensities. For the other observ- 
ables — Hi, IR, FUV, and Ha intensities — we use two- 
dimensional maps of intensity and determine the mean 
intensity for each tilted ring. 

Error bars on the Hi, IR, FUV, and Ha intensities 
show the lcr scatter within that tilted ring, capturing 
both statistical noise and deviations from axial symme- 
try. We estimate the lcr scatter for our CO measurements 
by integrating the CO cube over a velocity window that 
is adjusted for each line of sight such that it includes all 
channels of significant CO emission but at least all chan- 
nels with velocities within 25 km s _1 of the local mean 
Hi velocity. Note that our lcr values reflect the scatter in 
(integrated) intensities of individual lines of sight inside 
a ring. They should not be confused with the uncertainty 
in the determination of the mean intensity inside a ring 
which is typically much smaller. 

4. RESULTS 

In Figure H for NGC 628 and in the Appendix for the 
rest of the sample, we present stacked radial profiles of in- 
tegrated CO intensity along with profiles of Hi, infrared 
intensity at 24 fim and 70 /im, FUV, and Ha intensity. 
Following Equations (J3J) we combine Ha and 24 /im in- 
tensities to estimate the star formation rate surface den- 
sity SFR(Ha+24) and compare those to SFR(FUV+24) 
derived from FUV and 24 /xm intensities using Equa- 
tion (gj). 

For each galaxy there are two plots: The left panel 
shows Hi, H 2 (from CO), and SFR for both FUV+24/xm 
and Ha+24/xm. The right panel shows our SFR tracers 
- Ha, FUV, 24 /zm and 70 /an emission. We present 
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Fig. 4. — Azimuthally averaged intensity in 15" wide tilted rings in NGC 628: (left) Hi, H2, and SFR from FUV+24/im and Ho+24/^m; 
(right) 24 fim, 70 £im, FUV, and Ha intensity. Left-hand y— axes show observed intensities, the right-hand y— axis shows the surface 
densities of Hi, H2, and SFR projected to face-on values. The color of the CO points indicates the significance of the fit to the stacked 
spectrum: high quality in green, marginal quality in orange, and upper limits in red. The solid-dashed line shows an exponential fit to the 
radial CO profile with the scale length, ico, printed in the lower left corner. Error bars show the \u scatter inside each tilted ring. Note 
that we have scaled the Hi intensities (left y-axis) by a factor of 312.5 in order to match all profiles in units of surface density (see text). 



the profiles in both observed intensit^F 5 ! (left hand y- 
axis) and units of surface density (right hand y-axis of 
the left panel) — £h2, Shi, and SgpR. Note that we have 
projected (only) the surface densities of Hi, H2, and SFR 
to face-on values (i.e., we corrected for inclination). The 
observed surface brightnesses are not corrected for the 
effect of inclination. Such a correction will just move 
each galaxy up and down in lockstep and we find it more 
useful to report the observed values. The color of a point 
in the CO profile indicates the significance of the fit to 
the stacked spectrum: green for high significance, orange 
for marginal significance, and red for upper limits, these 
correspond to 3cr upper limits on the fitted intensity (see 
Section [33). 

4.1. CO and Star Formation 

With these azimuthally averaged data we are able to 
compare CO to tracers of recent star formation across a 
large range of H2-to-Hi ratios. In this subsection, we 
make empirical comparisons between measured intensi- 
ties, examine the relative roles of H2 and H 1 in the "star 
formation law," and investigate the origin of the scatter 
in these relations. 

4.1.1. Scaling Relations between CO and IR, FUV, and Ha 

Figures & [5] show scaling relations between observed 
intensities of CO (x-axis) and different tracers of recent 
star formation (y— axis). Figure [5] shows infrared inten- 
sities at 24 /jm (top panels) and 70 /im (bottom panels) 
and Figure [H] shows intensities of FUV (top panels) and 
Ha (bottom panels). The left hand panels show the rela- 
tions for all galaxies and all radii. The panels on the right 
hand side show only radii r > 0.5 r25- H2 and Hi make 
up roughly equal parts of the ISM near this radius (see 
the left panel of Figure [T3|) , so most of the points in the 
right hand panels are H i-dominated. A dotted vertical 
line shows an integrated CO intensity of 2.2 K km s _1 , 
which corresponds to Xh2 ~ 10 M pc~ 2 (assuming 
i = 45°), which is about the surface density at which 

15 In order to have Hi and H2 comparable in units of mass 
surface density we scale the observed H I intensity by a factor of 
312.5; the ratio of Equation {T} and J2j. 



H2 and Hi make up equal parts of the ISM. A dashed 
vertical line at Iqo = 0.3 K km s _1 (£h2 ~ 1 M pc~ 2 
for i = 45°) shows a conservative sensitivity limit for the 
whole sample. Typically our upper limits, which are not 
displayed in these plots, lie to the left of this line. They 
will be systematically higher at large radii when radial 
rings partially exceed the coverage of our CO maps. 

CO emission correlates tightly with IR emission at 
24 jim and 70 /im (Figure [5j (rank correlation coefficient 
fcorr = 0.9). The correlation extends over three orders 
of magnitude in CO and IR intensities and crosses the 
H i-to-H2 transition without substantial change in slope 
or normalization. Comparing the left panels (all radii) 
and the right panels (r > 0.5 ^5) does not reveal any 
significant radial dependence. 

CO emission exhibits a weaker correlation with FUV 
and Ha (Figure[6|) (r corr = 0.5 — 0.6) than with IR emis- 
sion, i.e., both the CO-FUV and CO-Ha relation show 
much larger scatter than the CO-IR relations. The CO- 
FUV relation displays a break between Hi- and H2- 
dominated regimes. The increased scatter and weaker 
correlation at least partially reflects the sensitivity of Ha 
and FUV emission to absorption by dust. In the inner 
(r < 0.5 ^25), more gas rich parts of galaxies dust re- 
processes most Ha and FUV emission into IR emission. 
In this regime CO and FUV are to first order uncorre- 
cted (r corr = 0.16). Outside ~0.5 r 2 5 the filling factors 
of dense gas and dust are lower and FUV is less affected 
by extinction. The correlation coefficient between CO 
and FUV or Ha is r corr = 0.45 — 0.55 in this regime, still 
not as high as for the whole galaxies because of the lim- 
ited dynamic range in intensity. In Section [4.1.2l we will 
see that the weaker correlation of CO with Ha and FUV 
is mostly due to galaxy-to-galaxy variations, possibly 
reflecting different star formation histories, dust abun- 
dances, geometries, and potential changes in the CO-to- 
H2 conversion factor. 

We use the ordinary least squares (OLS) bisector to fit 
power laws to each relation. In terms of sensitivity to 
a given amount of star formation, the SFR tracers are 
much more sensitive than our CO maps. To properly fit 
a relation between them we thus n eed to either caref ully 
incorporate upper limits (e.g., see iBlanc et al.1 120091 ) or 
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Fig. 5. — Observed IR intensities (j/-axis) as a function of integrated CO(2 — > 1) intensity (x-axis). Each point corresponds to a stacked 
average in a tilted ring 15" wide. Green and orange symbols indicate CO measurements of high or marginally significance; radial rings 
with only upper limits on the CO intensity (not shown here) are located exclusively to the left of the long-dashed line. The top panels 
shows the relation of 24 fim versus CO, the bottom panels shows 70 (im versus CO. The left panels show data for all radii, whereas the 
right panels show only data outside 0.5 r25, where the ISM is typically H I— dominated. The short-dashed vertical line indicates a typical 
CO intensity at which £h2 ~ Sjji; data to the left of this line will usually be H I— dominated. The diagonal dashed lines indicate lines of 
constant ratios for orientation. 



impose a matched sensitivity cut on the SFR tracer data. 
We take the latter approach, discarding data below 0.1, 
1.0, 10~ 3 , and 10~ 4 MJy sr" 1 at 24jim, 70 fim, FUV, 
and Ha intensities after an initial fitP°l. Graphically, this 
removes the flaring towards low intensity just above our 
sensitivity cut seen in Figure [5] & HI because we are not 
sensitive to a similar flaring towards low CO intensities. 

Table [5] reports these fits considering all regions with 
CO measurements of high significance and SFR tracers 
above the sensitivity cut. The main result is that CO 
emission is consistent with being linearly proportional 
to IR emission both at 24 /im and 70 fim, down to low 
surface brightness. The best fits relating CO with FUV 
or Ha emission are also consistent with a linear slope 

16 The result is relatively insensitive to the exact choice of sensi- 
tivity cut. Varying it by a factor of 2 affects the power law index by 
~0.05 for the CO-IR relations and ~0.15 for CO-FUV or CO-Ha. 



within the large uncertainties, but power laws are clearly 
an inadequate description of those data. 

4.1.2. Scatter in the Scaling Relations 

Each of the observed relations displays significant scat- 
ter: hi^m/Ico has lcr scatter of about 0.17 dex (~50%) 
and ho^ra/Ico scatters by 0.24 dex (~75%), whereas 
Ina/Ico scatters by 0.34 dex (a factor of 2.2) and 
Ifuv/Ico by 0.55 dex (a factor of 3.5). The origin of 
this scatter is of astrophysical interest. On small scales, 
this scatter arises from the evolution of individual star- 
forming regions, which vary dramatically in their ra- 
tios of CO-to-SFR tracers — leading to a breakdown of 
scaling relations when a resolution element c orresponds 
to an individual region (|Schruba et al.l f2010h . Here our 
azimuthal averaging washes out such small-scale varia- 
tions. Each point averages over many individual star- 
forming regions. However, we do have the ability to 
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TABLE 2 

Relation of IR Emission and SFR Tracers to CO 



SFR Tracer 


Rank Correlation 


Power Law Index 


24 fim 






... all data 


0.90±0.05 a 


1.0±0.1 b 


... r > 0.5 r 2 5 


0.87 ±0.10 


1.1 ± 0.1 


70 (im 






... all data 


0.87 ±0.05 


1.0 ±0.1 


... r > 0.5 r25 


0.84 ±0.10 


1.1 ± 0.2 


FUV 






... all data 


0.47 ±0.05 


0.8 ±0.6 


... r > 0.5 T25 


0.46 ±0.09 


1.1 ± 0.8 


Ha 






... all data 


0.63 ±0.05 


0.9 ±0.4 


... r > 0.5 T25 


0.55 ±0.10 


1.1 ± 0.6 



a We estimate the uncertainty in the rank correlation 
coefficient by taking the correlation coefficient derived 
from 1,000 random pairwise re-orderings of the data. 
b The slope quoted here is from the ordinary least 
squares bisector with the error estimated from the 
spread in fitting x vs. y and y vs. x. 



distinguish scatter within a galaxy from scatter among 
galaxies. 

To investigate the origin of the observed scatter, we 
remove galaxy-to-galaxy variations from the observed 
relation. We do so in two ways: First, we fit power-laws 
relating CO to IR, FUV, or Ha in each galaxy and then 
adjust all galaxies to have the same normalization. A 
priori we do not know if the galaxy-to-galaxy variations 
mainly affect the measurement of our gas tracer (x-axis) 
or the measurements of our SFR tracers (j/-axis) , there- 
fore, we repeat the exercise matching normalizations at 
a fixed value in y and then in x. We then compare the 
scatter among normalizations to the scatter about the 
re-normalized relation. 

We also carry out a more basic test, measuring how the 
scatter in the ratio of SFR tracer to CO emission varies 
both among and within galaxies. We compare the scatter 
in the median ratio of SFR tracer to CO, {hinm/Ico), 
among galaxies to the scatter in deviations from this me- 
dian ratio within galaxies, hi^m/Ico - (h^m/Ico)- 

For both approaches, we find that galaxy-to-galaxy 
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Fig. 7. — Relation between CO intensity and normalized FUV intensity (left panel) and normalized Ha intensity (right panel) for outer 
disks, r > 0.5 r25. For each galaxy we use power law fits to the data of high quality and normalize the relation by the FUV or Ha intensity, 
respectively, which the fit has at 1 K km s _1 . 



variations dominate the scatter in the observed relation. 
Scatter among galaxies in Figure [5] & [6] is ~2 times larger 
than the scatter within individual galaxies. The most 
striking cases are Ha and FUV emission in the outer 
parts of galaxies. Figure [7] shows that after normaliza- 
tion, the relation of CO with FUV and Ha emission for 
radii r > 0.5 r 2 5 becomes very strong (r colT = 0.9) and 
nearly linear (power law index of 0.90±0.05 for FUV and 
1.1±0.05 for Ha). Thus in the outer parts of galaxies the 
FUV-to-CO and the Ha-to-CO ratios vary dramatically 
among galaxies but are largely fixed inside each galaxy. 
The scatter appears driven at least in part by real sys- 
tematic variations in the ratio of CO-to-SFR tracers as 
a function of oth er galaxy parameters (see below and 
lYoung et al.lll996h . The case in the inner parts of galax- 
ies is more complex because of high dust attenuations 
resulting in non-linear CO-FUV and CO-Ha relations. 

4.1.3. H2 and Star Formation 

In Figure [8] we combine FUV and 24 /im intensities to 
estimate SsfRi which we plot as a function of Sh2- In 
Figure [9] we instead combine Ha and 24 /im to estimate 
Ssfr- In both Figures, the left hand panels show data 
for all radii, while the right hand panels show only rings 
with r > 0.5 T25, where the ISM is mostly Hi. As in 
Figure [5] & |H a vertical dotted line shows Eh2 ~ 10 
M Q pc -2 , a typical H 2 surface density where the ISM 
consists of equal part s H 1 and H 2 . 

In agreement with IBigiel et ahl ()2008l . 120 111) , we ob- 
serve an approximately linear scaling of SgFR and Efj2 in 
regions that are dominated by molecular gas (£h2 ^ 10 
M pc -2 ). For Figure|S]this agreement does not come as 
a surprise because we use many of the same data and a 
similar approach (FUV+24^m) to estimate recent SFR. 
The new results here are that (a) this trend continues 
without significant changes down to very low £h2, includ- 
ing regions strongly dominated by atomic gas, and (b) we 
find the same trend in Figure |H] using SFR(Ha+24/im). 

The rank correlation coefficient relating £h2 and Esfr 
for all data with at least a marginal CO measurement 
is r corr ~ 0.85. The scatter about a linear relation is 
0.25 dex (prior to any normalization). Both numbers 



resemble those derived for the CO-24 /im relation in 
Section 14.1.11 because 24 /im emission drives our hybrid 
SFR tracer over most of the area. In detail, the frac- 
tional contribution of 24 jxm to SFR varies with radius 
and choice of hybrid tracer. Generally speaking: (a) the 
larger the radius the larger the contribution of the un- 
obscured term, and (b) Ha contributes fractionally more 
than FUV to the hybrid. 

An OLS bisector fit yields a roughly linear slope and a 
molecular depletion time, Td C p = Sh2/Ssfr, of ~1.8 Gyr 
(including a factor 1.36 to account for heavy elements; 
see Table [3] for fit parameters) . This is sl i ghtly lower 
than Tdep = 2 G yr found by IBigiel et all (pOM) and 
iLerov et all ((2008) for a subset of the da ta analyzed 
here a nd Td cp = 2.35 Gyr recently found by IBigiel et al.l 
(|2011h for a sample that is similar to the one ana- 
lyzed here. We include (a) more starburst galaxies 
and (b) more low mass, low metallicity sp i ral ga lax- 
ies that IBigiel et al.l (|2008l ) and ILerov et all (|2008h ex- 
cluded from their CO analysis. Both dwarfs and star- 
bursts have shorter TH Rr > than large spirals (see below and 
Kennicuttl[T99llGao fc Solomoijl2004t ILerov et al.ll200l 
Daddi et al.ll2010f ). Moreover, our radial profiles weight 
these s mall galaxies more heavily than the pixel sampling 
used bv IBigiel et al.l (|2008l I2011D . 

As with the scaling relations between CO and tracers of 
recent star formation (Figure [5] & [6]) , the H 2 -SFR rela- 
tion exhibits significant scatter. We perform the same 
procedure to isolate galaxy-to-galaxy variations from 
scatter within galaxies and again find the scatter in the 
main relation (upper panels in Figure [S] & O dominated 
by galaxy-to-galaxy variations. We plot the relations 
after normalization in x and y in the middle and bottom 
panels of Figure [8] & [9] Once galaxy-to-galaxy scatter 
is removed, there is a remarkably tight, uniform linear 
relation linking molecular gas and star formation across 
almost three and a half orders of magnitude. 

What drives this galaxy-to-galaxy variation? In both 
the H 2 -SFR and observed intensity relations a large part 
of the scatter comes from a sub-population of less mas- 
sive, less metal-rich galaxies that exhibit high SFR- 
to-CO ratios. Figure [10] shows the metallicity depen- 



10 



Schruba et al. 



Observed Relation 



Observed Relation 



0.1 



0.01 



ict 



10" 



SFR(FUV + 24) 
all radii 



0.1 



0.01 



10" 



10" 



0.01 



10" 



10" 4 



4*5 
•OA 



c# 




0.1 1 10 100 1000 

Molecular Gas Surface Density [M Q pc~ 2 ] 

Relation normalized by y vs. x 



SFR(FUV + 24) 
all radii 




0.1 1 10 100 1000 

Molecular Gas Surface Density [M pc 2 ] 

Relation normalized by x vs. y 



SFR(FUV + 24) 
all radii 




0.1 1 10 100 1000 

Norm. Molecular Gas Surface Density [M G pc" 2 ] 



0.1 r 



0.01 



10 J r 



io- 4 L_ 

0.1 



SFR(FUV + 24) 
r > 0.5 r„ 



>> 0.1 - 



0.01 r 



10 J f 



10" 



0.01 r 



10" J f 



10" 



CP*' 



V 



1 10 100 1000 

Molecular Gas Surface Density [M pc~ 2 ] 

Relation normalized by y vs. x 



SFR(FUV + 24) 
r > 0.5 r«, 



VH 




0.1 1 10 100 1000 

Molecular Gas Surface Density [M© pc 2 ] 

Relation normalized by x vs. y 



SFR(FUV + 24) 
r > 0.5 r„. 



0.1 1 10 100 1000 

Norm. Molecular Gas Surface Density [M Q pc" 2 ] 
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Fig. 10. — Molecular gas depletion time, Td ep = ^H'i/^'.SKR: as 
function of metallicity. The depletion times shown here correspond 
to the ratio of the data shown in Figure[8]including regions with CO 
measurements of marginal significance. The median average T dep 
changes systematically with metallicity: for the low metallicity bin 
~0.8 Gyr, for the intermediate metallicity bin ~1.7 Gyr, and for 
the high metallicity bin ~2.4 Gyr. 

dence of the molecular depletion time, Td cp - Massive 
spiral galaxies with high metallicities have considerably 
longer (median averaged) depletion times: ~2.4 Gyr 
for systems with metallicities of 12 + logO/H > 8.75 
and -1.7 Gyr for 12 + logO/H = 8.65 - 8.75, while 
smaller galaxies with lower metallicities have systemati- 
cally shorter depletion times: —0.8 Gyr for systems with 
12 + logO/H < 8.653. We find a similar trend for the 
scaling of Td ep with the maximal rotation velocity, t!fl a t, 
that low mass systems with Vfl at < 140 km s _1 have 
Tdcp < 1 Gyr. These low mass, low metallicity systems 
are atomic dominated (Shi ^ Sh2) for most of their 
radii and show up prominently in the upper left panel of 
Figure El fc E as the dat a point s offset to shorter deple- 
tion times. iLerov et ail (|2008l ) and iBigiel et all (|2008ft 
labeled these galaxies "H i-dominated" and did not con- 
sider them in their H2-SFR analysis. 

4.1.4. SFR, HI, and H 2 

The question of which gas component — Hi, H2, or 
total gas — correlates best with recent star formation 
has received significant attention. Phrased this way, the 
question is not particularly well posed: the total gas sur- 
face density and the molecular gas fraction are closely 
related so that the different gas surface densities are not 
independent quantities. Therefore we do not necessar- 
ily expect a "best" correlation, only different functional 
forms. Still, it is illustrative to see how recent star for- 
mation relates to each gas tracer. Table [3] lists r corr and 
the power law index from an OLS bisector fit between 
each component and Figure [11] shows plots for each of 
the three gas phases and our two SFR prescriptions. The 
rank correlation coefficient and the power law fits are de- 
termined for regions where we have at least a marginal 
CO measurement. If we restrict our analysis to high sig- 

17 The absolute metallicity values that we quote should not be 
overemph asized, but the r e lative ordering of galaxies is fairly se- 
cure. See IMoustakas et al.l I pOlOh for more details. 



TABLE 3 

Relation of SFR to Different Gas Phases 



Gas Phase Rank Correlation Power Law Index 



SFR(FUV+24) vs. Gas Phase 


Hi 






... all data 


0.23 ± 0.05 a 




... Hi— dominated 


0.52 ± 0.06 




... H2-dominated 


0.41 ± 0.08 




H 2 






... all data 


0.88 ± 0.05 


1.0 ± 0.1 a 


... H 1— dominated 


0.75 ± 0.07 


1.3 ± 0.5 


... H2-dominated 


0.91 ± 0.07 


1.1 ± 0.1 


Hl+H 2 






... all data 


0.90 ± 0.05 


1.6 ± 0.2 


... H 1— dominated 


U. i ± U.UD 


n 1 n et 
Z.Z ± U.D 


... H2-dominated 


0.92 ±0.07 


1.2±0.1 


SFR(H 


ct+24) vs. Gas Phase 




Hi 






... all data 


0.21±0.05 a 




... H I-dominated 


0.54 ± 0.06 




... H2-dominated 


0.44 ± 0.08 




H 2 






... all data 


0.85 ± 0.05 


1.0±0.1 a 


... Hi— dominated 


0.70 ± 0.07 


1.5 ±0.7 


... H2-dominated 


0.91 ±0.08 


1.1 ±0.1 


H1+H2 






... all data 


0.90 ± 0.05 


1.7±0.2 


... H I-dominated 


0.79 ±0.06 


2.5 ±0.8 


... H2-dominated 


0.92 ± 0.08 


1.3 ±0.1 



a Estimation of uncertainties equal to Table [2] 



nificance CO measurements, we obtain the same results 
within the uncertainties. 

The table and figure show that recent star formation 
rate tracers rank-correlate approximately equally well 
with H2 and Hi±H2 both across all surface densities 
and separately in the Hi- and ^-dominated regimes. 
H 1 does not correlate significantly with star formation in 
the H2-dominated inner parts of galaxy disks, though 
the correlation between H I and recent star formation 
becomes stronger in the H I-dominated outer parts of 
galaxies (and furth er increases at even larger radii; see 
Bia iel et al.ll2010bh . 

The rank correlation is a non-parametric measure of 
how well the relative ordering of two data sets align. A 
high rank correlation coefficient implies a monotonic re- 
lationship but not a fixed functional form law. For total 
gas, the power law index relating gas and recent star 
formation depends fairly strongly on the subset of data 
used. If we focus on the regions where Shi > Sh2, the 
best-fit power law relating total gas and recent star for- 
mation (from FUV±24) has an index of 2.2 ±0.6. Where 
Ehi < Eh2, the index is much shallower, 1.2 ± 0.1. By 
contrast, the power law relating H2 to recent star for- 
mation varies less across regimes, from 1.3 ± 0.5 where 
Shi > S H 2 to 1.1 ± 0.1 where S ffl < S H 2- Both of these 
agree within the uncertainties with the fit to all data, 
which has slope 1.0 ± 0.1. The slight steepening of the 
relation in the outer disks (Shi > SH2) is driven by the 
interplay of two effects: the rather small dynamic range 
in gas surface densities and an increased dispersion in 
SFR-to-H2 ratios driven by the low mass galaxies; those 
have high SFR-to-H2 ratios and contribute mostly to 
the H i-dominated subset. The qualitative picture of a 
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Fig. 11. — SgpR (j/-axis) from FUV+24/xm (left) and Ha+24^tm (right) as a function of different gas phases: Hi alone (top), 
H1+H2 (micMZe), and H2 alone (bottom). Each point in these diagrams represents a radial average in a given galaxy. Regions that 
are H2-dominated are plotted with dark blue symbols, regions that are Hi— dominated in light red symbols. The bottom panels show only 
regions with at least marginal CO signal, while the top and middle panels show also regions where we determined only an upper limit on 
the molecular content. Whereas SFR is not correlated to Hi (in the inner parts of galaxy disks), it correlates with H2 and total gas. The 
scaling with H2 is uniform and linear for all regimes; the scaling with total gas exhibits a change in slope at the transition between Hi— 
and H2-dominated environments. 
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break in the total gas-SFR relation but a continuous 
H2-SFR relation remains unchanged if SFR(Ha+24) is 
considered, though there are small changes in the exact 
numbers (see the left and right panels in Figure 1111) . 

Our conclusions thus match those of Bigi el et al.l 
(2008): a single power law appears to be sufficient to 
relate £h2 and Esfr whereas the relationship between 
^hi+H2 and Esfr varies systematically depending on the 
subset of data used. Because we have a dataset that in- 
cludes significant CO measurements where Shi > £h~2 
we can extend these findings. First, total gas and H2 
are equally well rank-correlated with recent star forma- 
tion in all regimes and this correlation is always stronger 
than the correlation of recent star formation with Hi. 
Second, the H2-SFR scaling relation extends smoothly 
into the regime where £hi > £h2 whereas the total gas- 
SFR relation does not. Third, the result is independent 
of our two star formation rate tracers SFR(FUV+24) and 
SFR(Ha+24). 

4.1.5. Discussion of CO-SFR Scaling Relations 

Empirical Results: IR brightness at both 24 /im and 
70 fj,Ta correlates strongly with CO intensity over ^3 or- 
ders of magnitude. Across this range, there is a nearly 
fixed ratio of CO to IR emission. FUV and Ha emission 
show little or no correlation with CO intensity in the in- 
ner parts of galaxies, presumably due to extinction, but 
are found to correlate well with CO in the outer parts of 
galaxies after galaxy-to-galaxy scatter is removed. 

A result of these empirical scaling relations is that the 
ratio of recent star formation rate, traced either by com- 
bining FUV and 24 /j,m or Ha and 24 /im intensities, 
to molecular gas, traced by CO emission, does not vary 
strongly between the H i-dominated and ^-dominated 
ISM. This result is driven largely by the tight observed 
correlation between CO and 24 ^m emission. The tight 
relation between CO and 70 /jm emission suggests that 
the CO-IR relation actually holds for a larger range of 
mid-IR intensities. The tightening of the CO-FUV and 
CO-Ha relation in the outer part of galaxies (after re- 
moving galaxy-to-galaxy scatter) reinforce the idea of a 
linear relation extending to large radii. Thes e tight corre- 
lation s are consistent with the conclusion of iLerov et al.l 
( 2008) that there is only weak variation in the SFR per 
unit molecular gas mass with local environment. 

Galaxy-to-Galaxy Scatter: Each of the correlations we 
observe has significant internal scatter. Breaking this 
apart into scatter among galaxies and scatter within 
galaxies, we observe that in every case scatter among 
galaxies drives the overall scatter in the observed cor- 
relation. This "scatter" among galaxies is not ran- 
dom; less massive, less metal-rich galaxies exhibit a 
higher rati o of SFR tracer t o CO emission (see Fig- 
ure fTOl and lYoung et al.l[l99l . There are two straight- 
forward physical interpretations for this. The efficiency 
of star formation from H2 gas may be ge nuinely higher 
in the se sys tems, a view advoc ated by iGardan et al.1 
(|2007h and iGratier et all (|2010l ). Alternatively, CO 
emission may be depressed relative to the true amount 
of H 2 mass due to changes in the dust abundance. 
Low mass systems often have lower metallicities and 
correspon dingly less dust, which i s required to shield 
CO (e.g.. iMaionev fc Blackl H988t iBolatto et all H99l 
Idover fc Mac Lowll2010t iWolfire et al.H2010D . A precise 



calibration of the CO-to-H2 conversion factor as a func- 
tion of metallicity is still lacking, so it is not possible at 
present to robustly distinguish between these two scenar- 
ios. 

After removing these galaxy-to-galaxy variations we 
find a series of extraordinarily tight relationships between 
CO and tracers of recent star formation. The most strik- 
ing — and puzzling — example of this is the emergence of 
a tight correlation between CO and FUV emission in the 
outer parts of galaxies (r > 0.5 rag). This is puzzling be- 
cause one would expect the galaxy to be mostly causally 
disconnected over the timescales predominantly traced 
by FUV emission — 20 — 30 Myr compared to a dynam- 
ical (orbital) time of a few 100 Myr. Yet, somehow the 
differences between galaxies affect the CO-to-FUV ratio 
much more than the differences between the widely sepa- 
rated rings represented by our data points. Galaxy-wide 
variations in metallicity and dust abundance probably 
offer the best explanation for this. These may propa- 
gate into variations in the CO-to-H2 conversion factor 
and the average dust extinction. The latter may lead to 
scatter in estimates of SFR and both will affect Td cp - An 
alternative explanation is that external processes, which 
affect the whole galaxy, play a large role in setting the 
star formation rate on timescales traced by FUV emis- 
sion. 

These strong galaxy-to-galaxy variations partially ex- 
plain the unexpected lack of correl ation between CO 
emission and recent SFR observed bv iKennicuttl (pL998) . 
They averaged across whole galaxies and in doing so con- 
ceivably lost the strong internal relations that we observe 
but preserved the large galaxy-to-galaxy variations that 
offset internal relations. The result will be an apparent 
lack of correlation in galaxy-averaged data that obscures 
the strong internal relationship. Whether there is in fact 
a weaker relationship between H2 and SFR in galaxy- 
integrated measurements than inside galaxies depends 
on whether the suggested variations in star formation 
efficiency are real or a product of a varying CO-to-H2 
conversion factor. 

(The Lack of)A Molecular Star Formation Law: With 
improved sensitivity to CO emission we now clearly see 
nearly linear relations between CO and tracers of recent 
star formation rate spanning from the ^-dominated to 
Hi-dominated parts of galaxies. Note that the rela- 
tions will likely depart from these scaling relations if re- 
gions of high surface densities (Sh2 > 100 M@ pc~ 2 ) or 
starburst galaxies are considered (e .g., IKennicuttl 119981 : 
iDaddi et al.|[2010t iGenzel et al1l2010h . The lack of strong 
variations in the scaling between these two quantities in 
the "non-starburst regime" reinforces that molecular gas 
is the key prerequisite for star formation. A nearly lin- 
ear correlation over this whole range can also be restated 
as the absence of a strong relationship between the ra- 
tio Ssfr/Sh2 and Eh2- This implies that Eh2 averaged 
over a large area is not a key environmental quantity for 
star formation because it does not affect the rate of star 
formation per unit molecular gas. Apparently the global 
amount of H2 directly sets the global amount of star for- 
mation but the surface density of H2 does not affect how 
quickly molecular gas is converted to stars. 

By contrast, the host galaxy does appear to affect the 
ratio of star formation rate to at least CO intensity. This 
indicates important environmental variations but they 
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are not closely linked to surface density. In this sense, 
Figure [8] & [9] offer a counterargument against the idea of 
a star formation "law" in which gas surface density alone 
sets the star formation rate. Instead, over the disks of 
normal galaxies, we see star formation governed by two 
processes: (a) the formation of stars in molecular gas, 
which varies mildly from galaxy to galaxy but appears 
largely fixed inside a galaxy, and (b) the conversion of H I 
to H 2 , which does exhibit a strong dependence on envi- 
ronment inside a galaxy, including a strong dependence 
on surface density (Section 14. 2|) . In the second part of 
this paper we will look at this second process by measur- 
ing variations in the H2-H 1 balance as a function of gas 
surface density and radius. 

Systematic Effects: We have interpreted the ob- 
served scaling relations in terms of a relationship be- 
tween mole c ular gas and recent star formation rate. 
iLeroy et~ l. (2008) demonstrated that azimuthally aver- 
aged profiles of the FUV+24/im combination that we use 
here match those of several commonly used star forma- 
tion rate tracers with ^50% scatter. Here we have shown 
that using a combination of Ha+24/im to determine the 
SFR leads to indistinguishable results (see Leroy et al. 
2011, in prep, for more discussion). We now discuss 
several systematic effects that may affect the translation 
from observables to inferred quantities. 

The most serious worry is that the IR intensity, which 
drives the correlations, is acting as a tracer of dust abun- 
dance and not recent star formation. Gas and dust are 
observed to be well mixed in the ISM, so in the extreme, 
this would result in plotting gas against gas times some 
scaling factor (the dust-to-gas ratio). A more subtle ver- 
sion of the same concern is that CO emission is primarily 
a function of dust shielding against dissociating UV ra- 
diation. If there are large variations in the abundance of 
dust in the ISM then it may be likely that dust emission 
and CO emission emerge from the same regions because 
that is where CO can form and evade dissociation. 

A few considerations suggest that the 24 /an and the 
70 /zm emission are not primarily tracing dust abun- 
dance. First, over whole galaxies, monochromatic IR 
emission at 24 /urn and 70 /im does track the SFR (with 
some important var iations among types of galaxy; see 
iCalzetti et aT][20Tof ). Second, we observe a linear cor- 
relation with CO emission and not with overall gas col- 
umn, which one might expect for a dust tracer. Third, 
both the 24 /jm and 70 /jm bands are well towards the 
blue side of the peak o f the IR SEP for dust mixed with 
non star-forming gas (jBoulanger et al.|[l996l) and so are 
not likely to be direct tracers of the dust optical depth 
(mass). Still, a thorough investigation of the interplay 
between dust abundance, IR emission, and star forma- 
tion is needed to place SFR tracers involving IR emission 
on firmer physical footing. 

A less severe worry is that using dust and FUV emis- 
sion makes us sensitive to an old stellar component that 
might not have formed locally. Our targets are all ac- 
tively star-forming systems, so old he re means mainly old 
relati ve to Ha emission (r ~ 4 Myr; Mc Kee fc Williamsl 
119971 ). The appropriate timescale to use when relat- 
ing star formation and gas is ambiguous. When study- 
ing an individual region, it may be desirable to use a 
tracer with the shortest possible time sensitivity. Aver- 
aging over large parts of galaxies, one is implicitly trying 



to get at the equilibrium relation. Therefore a tracer 
with a somewhat longer timescale sensitivity may actu- 
ally be desirable. The typical 20 — 30 Myr timescale 
(jSalim et all 120071) over which most UV emission (and 
dust heating from B stars) occurs is well matched to cur- 
rent estimates for the l ifetimes of giant molecular clouds 
(|Kawamura et al. 2009) . This makes for a fairly symmet- 
ric measurement — with the spatial and time scales of 
the two axes matched — though one is comparing recent 
star formation with the material of future star-forming 
regions. 

The fact that SFRs derived from FUV+24/xm and 
Ha+24/zm are essentially indistinguishable indicates that 
the distinction between the time scale probed by Ha 
(~4 Myr) and FUV (20 - 30 Myr) is not important to 
this study, probably because of the large spatial scales 
considered by our azimuthal averages. 

There may also be systematic biases in our inferred 
Sh2 ■ We have already discussed the dependence of Vco 
on metallicity as a possible explanation for the high SFR- 
to-CO ratio observed in lower-mass galaxies. Vco cer- 
tainly depends on metallicity. Current best estimates 
imply a non-linear relationship, with Xg o sharply in- 
creas i ng below 12 + log in O/H ~ 8.2 - 8.4 jWoifire et al.l 
[20101 [Glover fc Mac Low! [20101: iLerov et all 120111 ). For 
our range of metallicities (~8.4 — 9.0) neither the es- 
timates of Vco nor the metallicity measurements are 
accurate enough that we feel comfortable applying a cor- 
rection to our data. Instead, under the assumption of a 
fixed SFR-to-H 2 ratio, the observed SFR-to-CO ratio 
can be utilized t o constrain the m e tallici ty dependence of 
the Vco factor. iKrumholz et al.l (|2011l ) adopt a version 
of this approach using literature data and show that the 
observed metallicity variation of the SFR-t o-CO ratio is 
broad ly consistent with Vco predicted by IWolfire et al.l 
(2010), though with large scatter. We therefore expect 
that Vco does affect our results, creating much of the ob- 
served offset to higher SFR for low-mass galaxies. Sub- 
sequent analysis, especially comparison to Herschel dust 
maps, will reveal if there are also important second-order 
effects at play within galaxies. 

Variations in the CO(2 — > l)/CO(l — > 0) line ratio cre- 
ate a second potential bias in Eh2 ■ We adopt a fixed ra- 
tio of 0.7 based on comparison to literature CO(l — > 0) 
data. This is somewhat higher tha n the observed rati o 
in the inner part of the Milky Way (Fixscn et al. 1999), 
~0.5, though th e uncertaintie s on t hat ratio are large. 
More important, iFixsen et al.l ([1999D suggest variations 
in the CO line ratios between the inner and outer Milky 
Way and there are well-established differences between 
normal disk and starburst galaxies. Although not im- 
mediately apparent from a co mparison of HERA CLES 
to literature CO(l -> 0) data (|Lerov et al.ll2009D . such 
variations could affect our derived Sh2 by as much as 
~50%. Rosolowsky et al. (2011, in prep.) will present a 
thorough investigation of how the line ratio varies with 
environment in HERACLES. 

4.2. Distribution of Molecular Gas 

The tight correlation between SFR tracers and CO 
emission across all regimes strongly reinforces the pri- 
mary importance of molecular gas to star formation. In 
this section, we therefore examine the distribution of 
molecular gas in galaxies. In the outer parts of spiral 
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Fig. 12. — (left) Normalized CO intensity (j/-axis) as a function of galactocentric radius (x— axis) in units of r25, the 25 th magnitude 
B— band isophote for a respective galaxy. The CO intensity is normalized so that an exponential fit to the profile has intensity 1 at 0.3 r25. 
Crosses mark the median and the la scatter among galaxies in a series of radial bins (including only galaxies with at least marginal data in 
that bin). The average decline can be parametrized by an exponential with scale length of 0.2 T25 with no clear evidence of a truncation or 
break, (right) Schematic diagram for an exponential disk with scale length 0.2 T25- We plot mean intensity, flux in each ring, and enclosed 
luminosity as functions of radius. Gray lines show r^o and rgo, the radii at which 50 and 90% of the flux are enclosed. 



galaxies where Shi > Sh 2 , the formation of molecu- 
lar gas from atomic gas appears to represent the bot- 
tleneck to star formation and the relative abundance 
of H 2 and Hi is key to setting the star f ormation rate 
(e.g-. lLerov et al.ll2008HBigiel et al.ll2010al) . We can ap- 
ply the large dynamic range in H2-to-H I ratios achieved 
by stacking to make improved measurements of how this 
key quantity varies across galaxies. 

4.2.1. Radial Distribution of CO Intensity 

Many previous studies have shown that azimuthally 
averaged CO emissi on decreases with increasing galacto- 
centri c radius (e.g.. lYoung fc Scovilldll991b | Young et al.l 
[1995L iRegan et all 120011 : iSchuster et all 120071 ). Whereas 
galaxy centers often exhibit deviations from the large 
scale trend, CO emission outside the centers declines ap - 
proximately uniformly with radius ([Young et al.l 1995). 
A first analysis of the HERACLES data revealed a char- 
acteristic exponential decline of CO emission in the inner 
parts of galaxy disks, with the scale length of CO emis- 
sion simila r to that of old sta rs and tracers of recent star 
formation ([Lerov et al.| [2009) . With the increased sensi- 
tivity from stacking and a larger sample, we can revisit 
this question and ask if this radial decline in CO intensity 
continues smoothly out to ~1 r 2 5. 

From exponential fits to the high-significance CO data 
of each galaxjE3 (solid-dashed lines in Figure [4] and Ap- 
pendix), excluding the galaxy centers (inner 30") we find 
a median exponential scale length, lco, of 0.21 r 2 5 with 
68% of all /co between 0.16 — 0.28 r 2 5. This value agrees 
well with typical scale l engths found in previous studies 
(e.g.. lYoung et al.lll995l ) and the individual galaxy scale 
lengths agree well with previous work on the HERACLES 

18 We exclude the following galaxies from the analysis because 
(a) their emission is compact compared to our beam: NGC 337, 
3049, 3077, 4625; (b) they are barely detected: NGC 4214, 4559; 
(c) their morphology is not well parametrized by an exponential: 
NGC 2798, 2976, 4725. 



sample (|Lerov et al.ll2009D . 

The normalizations of these fits reflect galaxy-to- 
galaxy variations in the total molecular gas content. In 
the left panel of Figure [12] we show all profiles aligned to 
a common normalization. We plot the radius in units of 
r25, the 25 th magnitude B-band isophote, and normalize 
each profile so that the exponential fits have intensity 1 
at r ga i = 0.3 T25- The figure thus shows the radial varia- 
tion of CO intensity across our sample, controlled for the 
overall CO luminosity and absolute size of each galaxy. 
Thick crosses mark the median CO intensity and the 
68th percentile range in bins 0.1 r 2 5 wide. The same ex- 
ponential decline seen in individual profiles is even more 
evident here, with Iqo — 0.20 ± 0.01 r 2 5 for the average 
of the sample. 

The left panel in Figure[T2lshows that the radial decline 
of the CO profiles observed previously for the inner part 
of galaxies extends without significant changes out to our 
last measured data points. In most galaxies there is no 
clear evidence for a sharp cutoff or a change in slope. As 
long as the normalized profiles are above the sensitivity 
limit the decline appears to continue (without significant 
deviation) with lco = 0.20 r 25 on average. This smooth 
exponential decline in CO intensity with increasing ra- 
dius suggests that the observed decline in star formation 
rate from the inner to outer parts of galaxy disks is driven 
by a continuous decrease in the supply of molecular gas, 
rather than a sharp threshold of some kind. 

There are several galaxies which deviate from this me- 
dian exponential trend. We already noted above that we 
do not fit exponential profiles to NGC 2798, 2976, 4725 
because their gas (Hi and CO) radial profiles are insuf- 
ficiently parametrized by exponentials. There are three 
galaxies (NGC 2146, 2903, 4569) that we do fit and de- 
termine small exponential scale lengths, Iqo ^0.1 r 2 5. 
NGC 2146 hosts an ongoing starburst and both NGC 
2903 and NGC 4569 have prominent bars that may be 
funneling molecular gas to their centers. For two other 
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Fig. 13. — The molecular-atomic ratio, i?H2 = Sh2/Shi as function of galactoccntric radius (left pan el), and total gas surfac e density, 
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galaxies (NGC 2841 and 4579) we determine large expo- 
nential scale lengths, ^co ^ 0-3 7*25. These galaxies are 
better described by a flat distribution (or even a central 
depression) and a cutoff at larger radii. 

The tight correspondence of the CO scale leng th, ^co 
to r25 has been noted before ([Young et al.lll995f h while 
other studies have found a close correspondence be- 
tween CO and near-infra red light (iRegan et al.l 120011 : 
iBlitz fc Rosolowskvl I2004D . In our sample there is a 
fairly good correspondence between r 2 5 and the near- 
infrared scale length, l^g, measured at 3.6 /urn with 
r 25 « 4.7 ± 0.8 l 3 . 6 (|Lerov et al.1 120081 a treatment of 
a larger sample suggests a slightly lower ratio of ~4.1 
with similar error bars). The near-infrared light should 
approximately trace the distribution of stellar mass, so 
that our measured scale length is very similar to that of 
the stellar mass. This tight coupling has been interpreted 
to indicate the importance of the s tellar potential well 
to co llecting star-forming material (jBlitz fc Rosolowskvl 
12006ft . Here we see this correspondence to continue into 
the regime where the molecular gas is not the dominant 
gas component, confirming that molecular gas formation 
is the bottleneck to star formation. 

The right panel in Figure [12] shows the distribution of 
enclosed luminosity, mean intensity, and flux at each ra- 
dius for an exponential disk with scale length of 0.2 r 2 §. 
The brightest individual ring for such a disk lies at 
0.2 T25 and half the flux is enclosed within r 50 ~ 0.3 r 25 . 
This value, r 50 , is fairly close to the radius at which 
Ehi ~ Sh2 in a typical disk gala xy (see the left panel 
of Figure 1131 and ILerov et al.ll2008ft . so that CO emission 
is about evenly split between the H2-dominated and H 1 - 
dominated parts of such a galaxy. Meanwhile, 90% of the 
flux lies within rgo ^0.8 r 2 § a valu e that is very similar 
to the threshold radius identified bv lMartin fc Kennicutt! 
(2001). We do not find evidence to support a true break 
at this radius, but as an "edge" to the star-forming disk, 
a 90% contour may have utility. 

4.2.2. The H 2 -to-HI Ratio 



In the outer parts of galaxy disks — and thus over most 
of the area in galaxy disks — we have Ehi ^ Eh~2 » i m ply" 
ing that star-forming H2 gas does not make up most of 
the interstellar medium. In this regime the relative abun- 
dance of H2 and H 1 is a key quantity to regulate the star 
formation rate. Observations over the last decade have 
revealed strong variations of the fraction of gas in the 
molecular phase as a function of galactocentric radius, 
total gas surface density, stellar sur face density, disk or- 
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two most basic of these trends, the H2-to-Hi ratio, 
i?H2 = Sjk/ShI) as a function of normalized galactocen- 
tric radius (left panel) and total gas surface density (right 
panel). We focus on i?H2 because it is more easily sepa- 
rated in discrete observables than the fraction of gas that 
is molecular, f m = SmA^m + S H2) = Rm/(1 + Rm)- 
This makes it easier to interpret uncertainties and sys- 
tematic effects like changes in the CO-to-Eb conversion 
factor. 

The left panel of Figure [13] shows i?H2 as function 
of galactocentric radius in units of r 2 5 for all data de- 
tected with high or marginal significance. For clar- 
ity, we do not plot i?H2 for regions where we deter- 
mined only upper limits in CO intensity; for the inner 
parts (r < 0.6 ^5) these upper limits are bounded by 
Rh2 ^ 0.3, whereas for outer parts (r > 0.6 ^5) the up- 
per limits in Eh2 are typically of comparable magnitude 
as measurements of Sri and upp er limits are bounde d 
by R m < 1 . In a greement with IWone fc Blitz! f2 0021 , 
iHever et al.l (|2004f ). lBigiel et al.l (|2008[ ). and ILerov et al.l 
( 2008) we find i?H2 to decline with increasing galactocen- 
tric radius, a variation that reflects the distinct radial 
profiles of atomic and molecular gas. However, radial 
variations alone do not explain the full range of observed 
molecular fractions because i?H2 can vary by up to two 
orders of magnitude at any given galactocentric radius. 

A significant part of the variations in i?H2 at a given 
galactocentric radius corresponds to systematic varia- 
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tions between galaxies. These are mainly caused by vari- 
ations in the absolute molecular gas content of a galaxy 
and can be removed by a normalization procedure sim- 
ilar to one applied in the left panel of Figure Q21 The 
result is similar to that seen for the radial profiles of CO: 
the relationship tightens and we can see that most of the 
decline of i?H2 inside a galaxy occurs radially. However, 
there is significantly more scatter remaining in i?H2 ver- 
sus radius than we observed for the normalized CO radial 
profiles, highlighting the importance of parameters other 
than a combination of radius and host galaxy to set i?H2- 

In addition to declining with increasing galactocentric 
radius, i?H2 increases as the total gas surface density, 
E gas = Shi+Sh2, increases. The more gas that is present 
along a line of sight, the larger the fraction of gas that 
is molecular. The right panel of Figure [T3] shows this 
result, plotting i?H2 as function of the total gas surface 
density, E gas . i?H2 increases with increasing E gas , with 
regions of high surface density, E gas > 20 M Q pc~ 2 , be- 
ing predominately molecular, i?H2 ^ 1< 

At high E gas , the right panel of Figure [13] is largely 
a way of visualizing the "saturation" of Ehi on large 
scales in galaxies. A number of authors have found 
that averaged over hundreds of p arsecs to kpc scales , 
S H i rarely exceeds ~10 M Q pc" 2 (|Wong fc Blitz! 12001 . 
This limit is clearly violat e d at high spatial reso- 
lution (e.g. iKim et all 119991 iStanimirovic et all 119991: 
Brinks k Shane 1984) and may vary among classes of 
galaxies. Several recent theore tical works have aimed 
at reproducing this behavior. IKrumholz et al.l ()2009D 
focused on shielding of H 2 i nside individu a l atom ic- 
molecular complexes, whereas lOstriker et al.l (|2010D ex- 
amined the interplay between large-scale thermal and 
dynamical equilibria. 

The right panel in Figure Q2] includes a lso the predicted 
solar- metallicity E gas -i?H2 relation from lKrumholz et al.l 
(|2009l their Equations 38 & 39). They model individual 
atomic-molecular complexes, however these complexes 
have a filling factor substantially less than 1 inside our 
beam. The appropriate E gas to input into their model is 
therefore the average surface densities of the complexes, 
E C omp, within our beam, which will be related to our 
observed surface density, E gas , by the filling factor, /, 
namely: E gas = / E comp . The model curve (blue short- 
dashed line) with a filling factor / = 1 is offset towards 
higher E gas (shifted right) or lower i?H2 (shifted down) 
compared to our data. To have the model curve intersect 
our data (red long-dashed line) requires a filling factor 
of / « 0.5 — 1, so that atomic-molecular complexes fill 
about half the areas in our beam. This is before any 
accounting for the presence of diffuse Hi not in star- 
forming atomic-molecular complexes and is assuming a 
fixed filling factor, both o f which are likely oversimplifi- 
cations. Nonetheless, the IKrumholz et al.l (|2009f ) curve 
does show an overall good correspondence to our data. 

As with the radius, the total gas surface density pre- 
dicts some of the broad behavior of i?H2 but knowing 
E gas does not uniquely specify the amount of molecular 
gas, particularly at low E gas . This are visible as the large 
scatter in i?H2 at low surface densities in the right panel 
of Figure Q2J The scatter reflects a dependence of i?H2 
on environmental factors other than gas surface density. 
Figure [14] shows i?H2 over the small range E gas = 7 — 12 
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Fig. 14. — The radial dependence of Rh2 for regions with com- 
parable total gas surface density, £hi+H2 = 7 — 12 Mq pc — 2 . The 
data are take n from the gray highlighted region in the left panel 
of Figure HHl Rn2 is n °t sufficiently determined by the total gas 
surface density but shows also a radial dependence. 

Mq pc -2 , i.e., the data from the gray highlighted region 
in Figure [13] We plot histograms for several radial bins 
which are clearly offset, indicating an additional radial 
dependence of i?H2- At small radii (r < 0.3 r 2 s) we ob- 
serve a large scatter in i?H2 for E gas = 7—12 M Q pc~ 2 
whereas at large radii (r > 0.6 r^s) gas with this surface 
density is always predominantly Hi. 

As was the case in SFR-H2 space, distinct populations 
of galaxies are responsible for some of the variations in 
the right panel of Figure [T3] Early type (Sab-Sb) spi- 
rals (e.g., NGC 2841, 3351, 3627, 4736) often show large 
molecular fractions, i?H2 ^ 0.5, but typically have low H 1 
and H2 surface densities, Ehi and Eh2 ~ 1 — 5 M Q pc~ 2 . 
By contrast, massive Sc galaxies (e.g., NGC 4254, 4321, 
5194, 6946) can have comparably high molecular frac- 
tions, i?H2 ~ 0.35 — 1, but have higher surface density 
Hi disks, Ehi ~ 5 — 10 M Q pc~ 2 , so that these fractions 
occur at higher E gas . A trend with metallicity is not 
immediately obvious in the data but these differences 
may reflect the more substantial stellar surface densi- 
ties found in the earlier-type galaxies. This increased 
stellar surface density results in a stronger gravitational 
field, which could lead to a higher midplane gas pressure 
(Blitz & Rosolowskv 2006) and a low fraction of diffuse 
Hi gas (|Ostriker et al.ll2010D . 

4.2.3. Discussion of Ru2 

Following several recent studies we observe strong sys- 
tematic variations in the H2-H1 balance across galaxies. 
Two of the strongest behaviors are an approximately 
exponential decrease in i?H2 with increasing galacto- 
centric radius and a steady increase in i?H2 with in- 
creasing gas surface density. Our improved sensitiv- 
ity shows these trends extending to low surface den- 
sities and our expanded sample makes clear that nei- 
ther of these basic parametrizations adequately cap- 
tures the entire range of i?H2 variations. The likely 
physical drivers for the scatter in i?H2 that we observe 
are metallicity and dust-to-gas ratio (e.g.. iLerov et al.l 
I2008t IKrumholz et al.l 120091 : iGnedin et al . 2009, Bolatto 
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Fig. 15. — The recent star formation rate per unit total gas, 
^SFr/^hi+H2 (the inverse of the total gas depletion time; y-axis), 
as function of the molecular— atomic gas ratio, i?H2 = ^Jb/^HI 
(x-axis). A fixed ratio of SFR-to-H2 is shown by the blue dotted 
line for a molecular gas depletion time of r dep ft: 1.4 Gyr and by 
the red dashed line for T^ ep f» 2 Gyr. In regions dominated by 
molecular gas (i?H2 > 1) the Ssfr/Shi+h2 ratio approximates 
a value corresponding to the constant molecular depletion time, 
T dcp- in regions of small molecular gas fraction (i?H2 < 1) the 
^SFr/Shi+H2 ratio decreases significantly indicating that the total 
gas does not scale uniformly with the recent star formation rate. 



et al. 2011, in prep.), the dissoc iating radiation field 
(e.g., iRobertson fc Kravtsovl 12008ft . variations in inter- 
stellar gas pressure and density (e g.. lEhnegreenl Il994t 
iWong fc Blit9 [20021: iBlitz fc Rosolowsky||200§ 7~and ex- 
ternal ^gerturjDation^tjh^t drive gas to higher densitie s 
(e.g., IBlitz fc Rosolowskvl [200l IBournaud etall [2010h . 
Each of these qu antities are observ ationally accessible in 
our sample fe.g.. ILerov et al.ll2008t ) and estimates of the 
relative roles of each process will be presented in Leroy 
et al. (2011, in prep.). 

The overall relationship of total gas and star forma- 
tion rate (middle panels of Figure lll[) can be reproduced 
by a roughly fixed ratio of SFR~to-H2 within galaxies 
(Section 14.1.41) a nd th e observed scaling of i?H2 with 
Ehi+H2 (Section I4.2.2|) . This is shown in Figure [T5l 
where we plot the SFR per unit total gas, Esfr/Ehi+H2 
(the inverse of the total gas depletion time) as function 
of the molecular-to-atomic ratio, i?H2 = Eh2/Ehi. At 
low surface densities, the i?H2 _ Effl+H2 relation regulates 
Esfr, at high surface densities where almost all of the 
gas is molecular the SFR-H2 scaling determines the ob- 
served ratio. This trend is well parametrized by a fixed 
SFR-to-H2 ratio with a molecular gas depletion time of 
Tdcp ~ 1.4 Gyr (blue dotted line) for all of our targets 
and Tdcp ~ 2 Gyr (red d ashed line) for big spirals (i.e. 
using the same sample as lBigiel et al.l 120081 ; ILerov et al.l 
2008) . The transition between these two regimes creates 
the curved shape seen in Figure [TT]&[l5j This offers more 
support for a modified versi on of the classical picture 
of a star formation threshold (|Skillmanl[T987l : iKennicuttl 
ll989UMartin fc Kennicuttll2001l ). in which dense, mostly 
molecular gas forms stars but the efficiency with which 
such (molecular) gas forms is a strong function of envi- 
ronment, decreasing steadily with decreasing gas surface 
density and increasing galactocentric radius. 



5. SUMMARY 

We combine HERACLES CO(2 -> 1) data with Hi 
velocity fields, mostly from THINGS, to make sensi- 
tive measurements of CO intensity across the disks of 
33 nearby star-forming galaxies. We stack CO spectra 
across many lines of sight by assuming that the mean 
H 1 and CO velocities are similar, an assumption that we 
verify in the inner parts of galaxies. This approach al- 
lows us to detect CO out to galactocentric radii ~1 r 2 5. 
Because we measure integrated CO intensities as low as 
0.3 K km s _1 (^1 M Q pc~ 2 , before any correction for in- 
clination) with high significance, we are able to robustly 
measure CO intensities in parts of galaxies where most 
of the ISM is atomic. 

Using this approach we compare the radially averaged 
intensities of FUV, Ha, IR, CO, and Hi emission across 
galaxy disks. We find an approximately linear relation 
between CO intensity and monochromatic IR intensity 
at both 24 /j,m and 70 /zm. For the first time, we show 
that these scaling relations continue smoothly from the 
H2-dominated to Hi-dominated ISM. Extinction causes 
FUV and Ha emission to display a more complex rela- 
tionship with CO, especially in the inner parts of galax- 
ies. In the outer parts of galaxy disks FUV and Ha emis- 
sion do correlate tightly with CO emission after galaxy- 
to-galaxy variations are removed. 

We use two calibration to estimate the recent star for- 
mation rate, FUV+24/im and Ha+24/im, which we com- 
pare to H2 derived from CO. We find an approximately 
linear relation between Esfr and Eh2 in the range of 
Eh2 ~ 1 — 100 M Q pc -2 with no notable variation be- 
tween the two SFR estimates. A number of recent studies 
(IBigiel et a I [20081 : ILerov et al.1 [200l [BlarjC^Tall [200l 
iBigiel et al.l 120 111 ) have also seen a roughly linear rela- 
tionship between Esfr and Eh2 and have argued that it 
implies that the surface density of H2 averaged over large 
scales does not strongly affect the efficiency with which 
molecular gas forms stars. 

We do find evidence for variations in the SFR-to-CO 
ratio among galaxies. Indeed, most of the scatter in 
the relations between CO and SFR tracers is driven by 
galaxy-to-galaxy variations. These variat ions are not 
rando m, but show the trend observed by lYoung et al.l 
(|1996ft that lower mass, lower metallicity galaxies have 
higher ratios of SFR-to-H2 than massive disk galaxies. 
It will take further study to determine whether these are 
real variations in the efficiency of star formation or re- 
flect changes in the CO-to-H 2 conversion factor due to 
lower metallicities in these systems. After removing these 
galaxy-to-galaxy variations the composite H2-SFR rela- 
tion is remarkably tight, reinforcing a close link between 
H2 and star formation inside galaxies. 

We compare the scaling between the surface densities 
of SFR and Hi, H 2 , and total gas (Hi+H 2 ). The re- 
lationship between SFR and total gas has roughly the 
same rank correlation coefficient as that between SFR 
and H2, but does not obey a single functional form. 
Where Ehi > Eh2 * ne relationship between Esfr and 
Ehi+h~2 is steep whereas where Shi < Eh2 the relation- 
ship is much flatter. Meanwhile, we observe a linear re- 
lationship between Eh2 and Esfr for the full range of 
Eh2 = 1 — 100 M© pc~ 2 . Ehi and Esfr are weakly cor- 
related and exhibit a strongly nonlinear relation, except 
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at very large radii. 

The unbroken extension of the £sfr _ ^H2 relation 
into the H i-dominated regime suggests a modified ver- 
sion ofJdie_dassfoal^^ 

old (ISkillmanlll987tlKennicuttl[T989l:IMartin fc Kennicutt] 
1200 ll) . in which stars form at fixed efficiency out of molec- 
ular gas, to first order independent of environment within 
a galaxy. The observed turn-over in the relation between 
SFR and total gas relates to the H 2 -to-H I ratio which is 
a strong function of environment. 

We therefore investigate the distribution of H2 traced 
by CO using our stacked data and compare it to the 
H 1. On large scales we observe CO to decrease exponen- 
tially with a remarkably uniform scale length of ~0.2 T25, 
again extending previous studies to lower surface densi- 
ties. We find the normalization of this exponential de- 
cline to vary significantly among galaxies. The H2-to- 
Hi ratio, traced by the ratio of CO-to-Hl intensities, 
also varies systematically across galaxies. It exhibits sig- 
nificant correlations with both galactocentric radius and 
total gas surface density and we present high-sensitivity 
measurements of both of these relationships. However, 
neither quantity is sufficient to uniquely predict the H2- 



to-H 1 ratio on its own. 
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APPENDIX 

Here we present the radial profiles of CO, Hi, FUV, Ha, and IR at 24 /xm and 70 /xm used to generate the plots in 
this paper. We average the data in 15" wide tilted rings. For the CO data, we stack the shifted spectra over this area 
and determine the integrated CO intensity from fitting line profiles to the stacked spectrum. For the Hi, FUV, Ha, and 
IR data we use two-dimensional maps of intensity (Section [3]) . We determine the la scatter from the 68 th -percentile 
from the data inside each ring. We plot these as error bars but note the distinction from the uncertainty in the mean. 

For each galaxy we present two plots: The left panel shows CO and H 1 both in units of observed intensities (K km s~ 1 , 
left-hand y-axis) and converted to mass surface densities (M© pc~ 2 , right-hand y-axis) of H2 and Hi. The color of 
the CO points indicates the significance with which we could determine the integrated CO intensities: green for high 
significance measurements, orange for measurements of marginal significance and red for 3<r upper limits. To have H2 
and Hi on the same mass surface density scale, we multiplied the observed 21-cm line intensities by a factor of 312.5 
(the ratio of Equations flU and ©). We also plot the star formation rate (SFR) surface density (M© yr _1 kpc~ 2 ) 
determined from Ha+24/im and FUV+24/im. Black solid-dashed lines show our exponential fit to the radial CO 
profile. We fit all high significa nce data excluding galaxy centers, d efined as the the inner 30", which often exhibit 
breaks from the overall profile dRegan et al.ll200H IHelfer et al.ll2003f) . The derived exponential scale lengths (in units 
of r25, the radius of the 25 th magnitude B-Band isophote), appear in the lower left corner. 

The right panel shows observed intensities (in MJy sr" 1 ) of our SFR tracers — Ha, FUV, 24 /zm and 70 /xm emission. 
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Atlas of radial profiles, see Figure[4]for details. 
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Atlas of radial profiles, see Figure[4]for details. 
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Fig. A. — Atlas of radial profiles, see Figure[4]for details. 
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Fig. A. — Atlas of radial profiles, see Figure[4]for details. 
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Fig. A. — Atlas of radial profiles, see Figure[4]for details. 
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